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Abstract 

We study the elastic response of a worm-like polymer chain with reversible kink-like 
structural defects. This is a generic model for (a) the double-stranded DNA with sharp 
bends induced by binding of certain proteins, and (b) effects of trans- gauche rotations in 
the backbone of the single-stranded DNA. The problem is solved both analytically and 
numerically by generalizing the well-known analogy to the Quantum Rotator. In the 
small stretching force regime, we find that the persistence length is renormalized due to 
the presence of the kinks. In the opposite regime, the response to the strong stretching is 
determined solely by the bare persistence length with exponential corrections due to the 
"ideal gas of kinks" . This high-force behavior changes significantly in the limit of high 
bending rigidity of the chain. In that case, the leading corrections to the mechanical 
response are likely to be due to the formation of multi-kink structures, such as kink 
pairs. 

PACS: 87.14.Gg, 87.15.La, 82.37.Rs, 87.15.Aa. 

1 Introduction 

Since the pioneering experiments of Smith et al. PP-[I], the single- molecule micromechanical 
studies have become one of the central themes in molecular biophysics. These techniques have 
contributed substantially into our understanding of the properties of various biomolecules and 
their interactions. The original stretching experiments with the double-stranded DNA (ds- 
DNA) were followed by the theoretical work of Marko and Siggia [5 a who demonstrated that 
the observed elastic response (i.e. the dependence of the ds-DNA end-to-end extension on 
the applied stretching force) is in remarkable agreement with the Worm-Like Chain (WLC) 
model of polymer elasticity. Within this model, the chain is described as a constant-length rod 
with finite bending rigidity subjected to the thermal fluctuations. Various modifications and 
refinements of the WLC model have been proposed over the past decade to include additional 
effects, such as torsional constrains [§], bond elasticity, and sequence disorder [HE]. Recently, 
a discrete version of the WLC has been suggested as a plausible model for the single-stranded 
DNA (ss-DNA) [Oj. In most cases, the modifications of the WLC model are within the domain 
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of linear elasticity, and the overall non-linearity of the observed response is associated with 
the entropic nature of the problem. 

In this paper, we discuss an intrinsically non-linear generalization of the WLC that takes 
into account effects of the localized structural defects on its elastic response. In particular, 
we study a generic model of the WLC with reversible kink-like singularities. This model 
is relevant to several systems of interest. First, the kink-like defects can be used to model 
the sharp bending of the ds-DNA backbone associated with binding of certain proteins |lUj - 
|21j . In addition, the protein- mediated looping of the DNA would normally result in a 
non-trivial mutual alignment of the ends of the loop, which can also be interpreted as a 
kink defect. Typically, such DNA-protein interactions are sequence-dependent, and there is 
a strong preference for the proteins to bind to the particular sites on the DNA. However, 
in a realistic situation the non-specific binding takes place as well, and it may affect the 
overall properties of the DNA in vivo. Since the stretching experiments can probe the elastic 
properties on the single-molecule scale, they may be used to extract the information about 
the parameters of such DNA-protein interactions. The theoretical study of this paper is 
important for the interpretation of the results of this type of future experiments. This 
problem was recently addressed in a general context of the DNA-deforming proteins [22] • I n 
that work, a discrete version of the WLC model was used to obtain the stretching curves 
with the numerical transfer matrix technique. In our paper, use of the continuous model 
allows us to obtain a variational analytical result and compare it against the numerics. In 
addition, our approach yields important insights into the underlying physics associated with 
the predicted elastic behavior. The only significant limitation of the continuous model is that 
it does not allow one to approach the (hypothetical) saturation regime where the DNA gets 
completely covered by the proteins. 

The other system where our model is relevant is the ss-DNA. While several models have 
already been proposed to describe its elastic properties, neither of them is sufficiently justified 
by the microscopic structure of the molecule, and all of these models break down at high 
enough stretching forces. Analysis of the backbone structure of the ss-DNA suggests that 
the traditional models such as the WLC or the Freely- Joint Chain (FJC) are unlikely to 
be adequate as the coarse-grained description of the chain. Indeed, on the atomic level the 
conformations of the ss-DNA are mainly associated with rotations of certain bonds. It is well 
known in the context of the general polymer physics that such rotations typically involve 
transitions between the discreet trans and gauche states, which are normally described with 
the Rotational Isomer model [2B] . On the other hand, there are also elastic modes associated 
with small deviations of the torsional angles from their local equilibrium values (i.e. from 
these trans and gauche states). The WLC would be a natural model for describing the 
long-wavelength elastic modes due to these small deformations. In order to account for both 
types of conformations, one needs to construct a hybrid description that would unify the two 
classical models of polymer physics. In fact, our problem of the WLC with reversible kinks 
has all the essential features of such a hybrid model. Indeed, the local trans conformation 
can be viewed as a "no-kink" state, while the gauche state can be considered as a kink since 
it corresponds to a turn in the overall direction of the backbone. 

Our model is different from a recent model by Wiggins, Phillips and P. Nelson [21], 
where the kinks are taken as freely-bending hinges instead of the fixed-angle singularities 
considered here. In most cases of the biological relevance the fixed-angle kinks appear as a 
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Figure 1: Kink geometry. 



more adequate description of the respective singularities than the free hinges. In particular, 
the gauche states are characterized by the specific angles between the bonds fixed by the 
chemistry of the polymer (ss-DNA). Similarly, protein binding and protein- mediated looping 
feature some characteristic angles determined by the nature of the proteins, although these 
angles may vary over some range of values (22] (we consider modification of our fixed-angle 
results for the case of these deformable kinks with thermal fluctuations at the very end of 
this work). Our results will be seen to depend on the kink angle quite substantially, and thus 
this extra parameter is important. A model with soft annealed kinks similar to [21] was also 
considered recently by Yan and Marko j2Hj in the context of DNA cyclization. 

In the following section we first describe our model for a worm-like chain with reversible 
kinks, and then outline a general approach to its solution by drawing an analogy to the 
Quantum Rotator. Subsequently, we solve the resulting evolution equation both analytically 
(by the variational method) and numerically (by direct integration), and compare the results. 
Finally, we discuss our results in both the weak-force and the strong-force limits as well as 
in the intermediate regime, and draw physical conclusions from our findings. 

2 Theory 

Model. We consider a worm-like polymer chain with rod-like bending elasticity. The 
persistence length of the unperturbed chain is denoted by l p , and the total length of the 
chain L is much greater than the persistence length throughout this paper. We presume that 
kink-like structural defects can exist anywhere along the chain. These defects are reversible 
and can appear and disappear spontaneously, with certain free-energy penalty associated 
with each of them. We should emphasize that we neglect the sequence dependence of kink 
probability, thus limiting ourselves to the case of non-specific DNA-protein binding only. 
Each kink is characterized by the (fixed) opening angle K (Fig. QJ; this value is the same 
for all the kinks and is an external parameter of the problem. The presence of each kink 
costs finite amount of the free energy e (at zero stretching force), which is also assumed to 
be the same for each kink. The structural defects are local, i.e. they are characterized by a 
microscopic length scale Iq much smaller than the persistence length l p and can be considered 
as point-like for most practical purposes. We also neglect any direct interaction of the kinks 
separated by more than l . Thus, the three length scales of the problem are related by 
L > l p > l . 

An important quantity is the average line density of the kinks k in the absence of the ap- 
plied stretching force. This quantity is defined as the Boltzmann probability of the existence 
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of a kink per characteristic length scale of a kink site 1$: 



n= 1 exp(-e/A;T) 
Iq 1 + exp(— e/kT) 

At each site of length Z 0; there might be either a kink [factor exp(— e/kT)] or no kink [factor 
1], with the probability of having a kink at an arbitrary site (per unit length) given by the 
expression above. Note that in practice we assume that e 3> kT, and therefore 

k « r exp(-e/A;T) < (2) 
to *o 

i.e. kinks are rare and far apart, and only few possible sites are occupied by the kinks. 

The chain is stretched by applying force F in the z direction of the Euclidian coordinates. 
The effective energy of a chain segment between two adjacent kinks is given by the sum of 
the bending energy and the coupling to the stretching force: 



E f S i+l 

~kT = 




ds. (3) 



Here s is the coordinate along the chain, t(s) is the unit tangent vector of the chain at s, 
and kink number i is located at coordinate s, [so that L = J2i( s i+i ~ s i))- This expression 
was used by Marko and Siggia [5 j and in some earlier works on worm-like chain elasticity. 
This time, however, it does not apply to the entire chain; instead, it applies only to the 
chain segments between kinks. At each kink (i.e. at s$) the tangent vector experiences an 
abrupt change of its orientation and hops to a direction n — K away from the preceding one. 
Mathematically, this can be written constraint 

t( Si - 0) ■ i( Si + 0) = cos(tt -K) = - cos K. (4) 



Analogy to the Quantum Rotator, and the free energy. Let tp(i, s) be the ori- 
entational distribution function for the chain ends. Then the probability distribution for the 
tangent vector t at coordinate s for inner points of the chain is simply \ip(t, s)\ 2 . Function ip 
satisfies a Schrodinger-like equation for evolution along the chain. This can be demonstrated 
by the path-integral technique I2Z1 where path integration is conducted over all possible 
chain conformations, thus accounting, among other things, for the entropy of the chain. A 
more illustrative (though less rigorous) way to obtain the same equation also exists [2*%] . 
Here, we will only notice that each term in the effective energy (jSJ) has a counterpart in this 
Schrodinger-like equation, in the same way as each term of the mechanical energy in classical 
mechanics has a counterpart in the quantum-mechanical Hamiltonian. Thus, an unperturbed 
chain (without any kinks) obeys the evolution equation d s ip(t,s) = Hoi/)(t,s), where Hq is 
the effective Hamiltonian 

The term with the Laplacian in t arises from the term with the s-derivative in the effective 
energy (J3J), while coupling to the external force remains the same [up to the sign change] (see 
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Refs. [5J l2Zj for details). The Hamiltonian H is the Hamiltonian of the Quantum Rotator 
in quantum mechanics. 

Now we need to add a term to this Hamiltonian responsible for the non-local constraint (j3J). 
This teleportation-like term changes the direction of t abruptly, with the tangential vector 
t(s) hopping instantly by a finite angle 7r — K at the location of a kink. Such behavior is 
described by a delta-functional kernel in the Hamiltonian: 

r d 2 i' 

VV>(t,s)oc / — <?(t-t' + cosKW>(t / ,s), (6) 

J Z7T 

where the full Hamiltonian is now H = Hq + V and the 1/2tt factor provides the proper 
normalization of the kernel. The number of such non-local hops of the tangential vector (i.e. 
the number of kinks) yields the proportionality constant in the expression above. Thus, the 
missing pre-factor is simply the average line density of the kinks k, and the full Hamiltonian 
is 

Hip{i, s) = 2]-A { V>(t, s) + ^z- tV>(t, s) + kJ ^5{t ■ t' + cosK)iP{i', s). (7) 
The structure of the evolution equation 

^ = a*M, 00 

suggests the standard procedure for solving this kind of quantum-mechanical-like problems: 
expansion in eigenf unctions of the time-independent (i.e. s-independent) Schrodinger equa- 
tion 

-/^(t)=#V(t), (9) 

where ^(t) and \i are the eigenfunctions and the eigenvalues, respectively. Since each of the 
terms in such an expansion depends on the exponent of (the negative of) the corresponding 
eigenvalue, it is easy to conclude that for sufficiently long chains most terms in the expansion 
die off very fast, and only the smallest-eigenvalue term governs the long-chain behavior. 
Thus, we are interested in the lowest eigenvalue \i of the equation 0. The free energy of the 
chain is then simply related to this lowest eigenvalue by 

w = L » < 10 > 

(see Refs. 0EZ| for details). 

Variational solution. Since the lowest eigenvalue is being sought, we use the varia- 
tional method for finding the analytical solution to eigenproblem Q. Application of the force 
makes the z direction highly preferable for the tangential vector, and therefore the distribution 
function ^(t) must be strongly-peaked at t||z. Conventionally, we choose ^(t) oc exp(a>z-t/2) 
as a trial function, or 

I UJ ( uz ■ t\ 
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upon proper normalization (/ ip 2 (i) d 2 i = 1). Here uj is the variational parameter. The lowest 
eigenvalue can then be estimated as 



maxiJ, 



'121 



where 

h = J ^(i)Hip(i)d 2 i. 

Evaluating the last expression with trial function (fTT]) yields 



mm 




cothu; 



K 



sinh(au;) 
a sinh uj 



(13) 



(14) 



where we introduced a = sin(K/2). Minimization with respect to uj gives the implicit 
dependence of uj on the applied force or, being inverted, the explicit dependence of the force 



on uj: 



F uj 2 cosh uj sinh uj — uj 
kT 4L sinh 2 uj — uj 2 



+ KM* 



,Ot 



_1 sinh (aw) cosher — cosh (aw) sinh a; 



sinh uj 



uj< 



(15) 



Note that here and everywhere below uj is the particular value of the variational parameter 
that minimizes the free energy (instead of the generic variable implied until this point). 
Knowledge of the free energy as a function of force ()14j) also allows one to determine the 
extension of the chain: 



z 1 dT 



L LdF d(F/kT) 



coth uj — 



1 



to 



(16) 



Thus, equations (|T5|) and (|T6|) provide the parametric dependence of the force on the exten- 
sion, with uj playing the role of the parameter. This extension curve is shown in Fig. 121 for 
kink angles of 135° and 45° and several values of k. Our results reduce to those of Marko and 
Siggia [S] when no kinks are present, i.e. when k — 0. The new physics introduced by the 
kinks is represented by the term proportional to k in equation ()15j) . The ratio of the first and 
the second terms in this equation is determined by the dimensionless parameter nl p , which 
has the physical meaning of the average number of kinks per the persistence length in the 
absence of force. The larger this parameter is, the more significant the kink contribution is. 



Weak forces and small extensions. The asymptotics of the above results in the limit 
of small uj: 

and 



F 
kT 



2T p + < 1 + a2 



uj + O (uj 



allow one to obtain the explicit dependence of the extension on the force in the limit of weak 
forces and small extensions: 



L 31 + 2kL(1 — a 2 ) kT' 



(19) 
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Figure 2: Extension curves: force vs. extension. Analytical results of the variational solution 
(a) for K = 135°, and (b) for K = 45°. Curves for several values of the kink density kI p are 
shown on each plot. 
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This expression can be compared to the weak-force result for a bare (kink-free) persistence 
chain z/L = (2/3)l p (F/kT). Clearly, the elastic response of the chain with kinks is charac- 
terized by the renormalized persistent length 

lpeS = 1 + 2kZ p (1 - a 2 ) = 1 + 4,(1 + cos if)' ^ 20) 

This quantity reduces to the bare persistence length l p in the absence of kinks (k — 0). It 
is also apparent that l pe g < l p , and therefore a worm-like chain with kinks is harder than 
a worm-like chain without them. Similar renormalization of the persistence length was also 
observed in recent numerical results by Yan and Marco [22] • It may appear that analogous 
renormalization was obtained in Refs. [3 |H| as well, however, these works consider quenched 
sequence disorder, and thus the physical origin of their renormalization is very different. 

For high bending rigidity (kI p 3> 1), the renormalized persistence length becomes l pc s = 
l/[rc(l + cos if)]. This result deserves some discussion as it turns out to be closely related 
to the Flory result for trans-gauche rotational isomers [2H|- If a gauche conformation is 
considered as a kink, then it is straightforward to identify the geometrical relations 

K V3 . 6 
cos — = — sm- (21) 
2 2 2 K J 

and Iq = 6cos(#/2), where 9 is the angle characterizing the trans zigzag and b is the length 
of the zigzag segment (Fig. EJ). From these results, one can obtain the mean square of the 
distance between the chain ends: 



■W A7 2_ 2 2 + ex P {jr) ■ otf 
k ° 3sin 2 (f) 



(R 2 ) = 2l pcS L = 2^Nl 2 = 2 - , ,}^ N9 cos 2 v -, (22) 



which is identical to the result for the Flory model 

1 — cos v 3 

in the limit of rare kinks (e ^> kT). Note that comparison to the Flory result for the free 
rotational model (instead of the trans-gauche rotational one) would be impossible due to the 
lack of the equivalent angle and absence of a relation similar to Eq. (|2ip. In the case of the 
free rotations, the effective kink angle can adopt a continuous spectrum of values between 
zero and some positive value corresponding to the rotation by 180°. 

Strong forces and large extensions. In the opposite limit of large ui our results (|15|) 
and (J16|) also allow for simple asymptotic expansions: 

z 1 . . 

7 = 1-- (24) 

and 

¥i = w p + exp [-( ~ a)uj] • ( } 
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Figure 3: Geometrical parameters of the Flory model for trans-gauche rotational isomers [23*] . 



Thus, an explicit dependence for the elastic response is again available: 



z 
L 



\ 



kT 
FL 



2 p a 



a 



exp 



-2(1 




(26) 



To the main order, this result reproduces the earlier result of Marko and Siggia jH] for the one- 
over-the-square-root dependence on the stretching force. Hence, the response of a worm-like 
chain with kinks to the strong stretching is dominated mostly by the bare worm-like-chain 
elasticity In addition, kinks contribute an exponential correction to that main-order result, 
which can be interpreted as due to the "ideal gas of kinks". Indeed, the second term in 
equation (jl4j) is the average line density of the kinks in the presence of the stretching force, 
and it is equivalent to the concentration of the molecules in the ideal gas. Thus, the second 
term in equation ()25|) derived from the free energy (jl4j) is analogous to the relation p = nkT 
for the pressure of the ideal gas. 



High rigidity, small kink angle, numerical results. Probably the most interesting 
feature of Fig. El is related to the middle portion of the high-rigidity extension curves, where 
an instability region with dz/dF < is present. This instability region exists when parameter 
kI p is greater than the "critical" value (kI p ) c , i.e. when the chains are of high rigidity. The 



"critical" value is determined by the conditions F'(z) = F" 



and is given by 



(kI p ) 



exp uj c 



sinh aui c 



(4u; c (1 + a 2 ) - 8) + cosh au c (8 - 8u; c ) 



(27) 



where o> r is the solution to 



tanh aoo 



cu(3 + a 2 ) -6 



a 



lo(1 + 3a 2 ) - (3 + 3a 2 ) 



(28) 



and is large for all kink angles (the lowest uj c is around 5 for small kink angles). The last 
equation allows for simple analytical solutions in the limits of small kink angles (au 1) 
and large kink angles (au ^> 1). Thus, the analytical asymptotics of (k/ p ) c in these two 
limits are readily available: 



{kI p ) 



exp (3 + v^) / 
8(1 + V3) V 



+ 3a 



,2 + y/3 N 
3 + y/3, 



{au < 1) 



(29) 
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Figure 4: "Critical" value of parameter kI p as a function of the kink angle K. Both the 
numerical solution and the analytical asymptotics are shown. The instability region is above 
the solid curve. 



and 

. , . aexp3 . , . 

Wc = oh , (au » 1) . (30) 

Z[l — a) 

These asymptotics as well as the exact numerical solution for {kI p ) c are shown in Fig. 0] As 
can be observed in this figure, the "critical" value is much lower for sharp kinks, and hence 
for smaller kink angles the instability region is present for much lower values of kI p and lies 
in the range of the measurable stretching forces in Fig. El 

This instability behavior is similar to the behavior of the Van-der-Vaals curves for the 
ordinary liquid- gas system and hence suggests the co-existence of two "phases" . Of course, no 
true phases can exist in a one- dimensional system like the DNA molecule under consideration, 
however, for the worm-like chain with kinks it may be appropriate to think of "kink-free" 
and "kink-rich" portions of the chain instead of the true phases. 

In order to check this "two-phase" hypothesis, we conduct a direct numerical solution of 
the evolution equation (JHJ) with the Hamiltonian (J7J). In order to avoid the delta-functional 
singularity inside the integration in the last term, the analytical integration of the delta- 
function is carried out, and we operate with a smooth-kernel Hamiltonian 



Hip{p,s) = — 



l _ ^ d 2 ^(p, s) _ 2p dip(p, 



dp 2 ' dp 



F 

+ —pi){p,s) + 



K r~P cos K+^l-p 2 sinK lb(p' s) dp' , , 

+ - / . , 1 ; (31) 

7T J-pcosK-^l-p^inK J (I -pi) (1 -p'2) - (pp> + COS K) 2 



(where p — z • t and p' — z • t') instead of the original Hamiltonian (|7J). Then we choose 
an arbitrary initial distribution function ip(t, 0) and let it evolve by evaluating the right- 
hand side of Eq. (jHJ) on each step and setting the s-derivative (the left-hand side) equal to 
the result. This procedure is stopped when the relative change of the relative extension 
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z/L — J z • iijj 2 (i } s) d 2 i evaluated at each step becomes very small and hence s) itself 
becomes the stationary distribution function. Such an approach yields more precise results 
for all the quantities since it does not make any a priory assumptions about the shape of the 
distribution function. 

The results of this numerical solution are shown in Fig. El for kink angles of 135°, 90°, 
and 45° and several values of k. It is immediately apparent that the agreement between 
the variational analytical and the exact numerical curves is good for both small and large 
forces (and extensions) for all values of parameters, and hence the analytical results of the 
preceding two sections are accurate in the respective limits. However, the intermediate 
regime for high-rigidity chains and small kink angles differs substantially from the analytical 
parametric dependence of equations (fTK|) and QlfiJ). What causes this difference, and how can 
the numerical results be understood? 

The answer to both these questions comes from the numerical results for the distribution 
function ip(t,s) (Fig. [HJ) - The distribution profile for large kink angles possesses the shape 
assumed in the variational method: a peak at t along z with the exponential fall-off away from 
it. However, for rigid chains and small kink angles, we observe existence of a secondary peak in 
addition to the main peak. Thus, the discrepancy between the variational and the numerical 
results can be explained by inadequacy of the trial function in the variational method, which 
ignores the possibility of the secondary maximum. The location of the secondary peak 
corresponds approximately to the angle tt — K away from the direction of the applied force 
(z), i.e. there is an unusually high portion of chain segments at angle tt — K to the force. 
This can be interpreted as presence of a small, but substantial fraction of kink pairs (Fig. EJ) • 
Indeed, if the two outer chain sections of a kink pair are aligned precisely along z, then the 
inner section is at angle tt — K to that direction, and if the kinks in each pair are close, 
then the fraction of the chain segments at angle tt — K to the force is small compared to the 
fraction of the chain segments aligned with the force, but high compared to the fraction of 
the other non-aligned segments [in full agreement with Fig. Efb)]. When the two kinks are in 
a pair with the outer sections aligned with the force and the inner section at angle tt — K to 
the force (Fig. EJ), much less bending is necessary compared to the case of a single kink with 
both outer sections aligned with the force (Fig. Thus, condensation of kinks into pairs 
is favorable both because the fraction of the "non-aligned with the force" chain segments 
is small and because the energetically costly bending is not required. We should emphasize 
that although kink-pairing may be a good way of thinking of this system, we do not have a 
simple yet rigorous theoretical model accounting for all the properties of such a "gas of kink 
pairs" . 

Soft kinks. So far, we considered only hard kinks, i.e. kinks characterized by the fixed 
opening angle K. Each kink could appear and disappear; however, the possibility of its elastic 
deformation due to the thermal fluctuations was ignored. In this section, we demonstrate how 
our results are modified when kinks are of finite rigidity and their opening angle is allowed to 
deviate from the equilibrium value K. We limit our discussion to small fluctuations around 
K, which corresponds to high rigidity of each kink. In the opposite limit of large deviations, 
the kinks are almost free and the model of Wiggins et al. (21] should be applied. 

The effective energy of the chain can be written as a sum of the contributions (jSJ) of all 
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Figure 5: Extension curves: force vs. extension. Numerical results of the direct solution of 
the evolution equation (a) for K = 135°, (b) for K = 90°, and (c) for K = 45°. Curves for 
several values of the kink density kI p are shown on each plot (same values as in Fig- EJ) - Solid 



lines represent the numerical results, 
analytical results. 



while the dashed lines reproduce the corresponding 
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Figure 6: Numerical results for the distribution function ip(p, s) (where p — z • t) in the limit 
of large s for the kink angles of (a) 135° and (b) 45°. The two graphs are for the comparable 
stretching force and for the respective values of the kink density [\^10 times the "critical" 
density for each kink angle, i.e. for the second from the right curve in Figs. Efa) and^c)]. 




Figure 7: A pair of kinks. 
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segments between kinks and the quadratic contributions of all kinks: 



E. 



soft 



kT 



E 




kT 



-z • t 



ds + 



7 



arccos 



(-t( Si -0)-t( Si + 0)) -K 



(32) 

where the summation is over kinks i. The last spring-like term allows for the thermal fluctu- 
ations of the opening angle and substitutes for the delta-functional constraint (@J); a similar 
term was used in the numerical study by Yan and Marko [22J. Kink stiffness parameter 
7 is assumed to be large, so that the fluctuations are small. The V'-term in the effective 
Hamiltonian is then 



KoftVKM) = K 



-ii')-K\ 



iP(i',s)d 2 i' fe 



.X(K'-K) 2 



V?p(i,s)d cos K' 



fe 



-tt' 



-K 



d 2 l' 



j e -m'-K) 2 dcosK' 



(33) 



where in the last equality V is the old delta-functional term (JHJ) (multiplied by k) considered 
as a function of the integration variable K' instead of the equilibrium value K. 

Thus, the variational solution for the free energy of the chain with soft kinks reduces to 
a single integration of the last term of the variational result ()14|) for the free energy of the 
chain with hard kinks. This integration can be performed analytically in the limit of small 
thermal fluctuations (7 ^> uj/a) and leads to the following main order result: 
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or, in terms of the force, 
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(35) 

In the expressions above we omitted all the terms of the order of uj/'j or less and retained only 
the term of the order of uj 2 /^, which is not necessarily small when uj is large. Thus, the main 
effect is simply renormalization of the kink density k by a factor of exp [o; 2 (l — a 2 )/ (87)]. 
Physically, when kinks are allowed to relax by deforming the opening angle, the number of 
kinks increases, making the chain harder to extend. However, the absolute value of the effect 
is not substantial, and for all reasonably small thermal fluctuations of the kink angle (large 
7) the analytical extension curves get only slightly distorted for large extensions, without 
major qualitative changes. For very soft kinks (small 7) the results of Wiggins et al. |21j 
should be employed. 



3 Conclusions 

In conclusion, we studied the generic model of the semi-flexible polymer chain with reversible 
kinks. It can be viewed as a hybrid of the two classical descriptions of the polymer elasticity: 
the WLC and the rotational-isomer-states models. Therefore, the proposed theory should be 
applicable to the ss-DNA whose conformations may involve both the discrete trans-gauche 
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rotations of the chemical bonds and small deformations giving rise to the continuous elasticity. 
Another important class of systems where our model is applicable is the ds-DNA with sharp 
protein-induced bends. 

In the limit of weak stretching forces, the elastic response of the kinked DNA chain is 
characterized by a renormalized persistence length, which is smaller than the bare persistence 
length. This conclusion is consistent with the observation made in the earlier numerical work 
on the problem [22J. We obtained the analytical expression for the renormalized persistence 
length and showed that the classical results for both the pure WLC and the rotational 
isomer models can be recovered exactly as limits of our expression. Furthermore, by using 
the variational approach, we calculated the complete non-linear response of the chain to the 
stretching. This result is in excellent agreement with the direct numerical solution over a 
substantial range of the model parameters. 

In the limit of strong stretching forces, we recover the pure worm-like-chain behavior with 
exponential corrections due to the "ideal gas" of kinks. The variational theory breaks down 
in the regimes of high chain rigidity and small kink angles. In this case, the analytical curves 
have signatures of instability similar to those of the Van der Waals gas. By analyzing our 
numerical results, we conclude that this behavior corresponds to the creation of multi-kink 
objects, e.g. kink pairs. 

For soft kinks, where the opening angle can fluctuate and thus relax the overall energy of 
the chain, we found that the number of kinks increases compared to the hard-kink case, and 
the DNA gets harder to extend. 

The major limitation of our model is that we have neglected the sequence-specific effects 
by assuming that the kink energy is constant along the chain. While this may be a reasonable 
first approximation to both problems of the ss-DNA elasticity and the non-specific protein- 
DNA binding, a significant future work is needed in order to include the effects of the sequence 
disorder. 

Several experimental tests of our results can be suggested. In the case of the ss-DNA, 
in order to probe the "pure" elastic response of the chain one needs to exclude the effects 
of base-pairing and the electrostatic interactions. In the existing experiments, these effects 
are not suppressed, and thus direct comparison is not possible at the moment. (Note that 
the electrostatic and the base-pairing contributions were simulated in Ref. [2E] instead of 
being excluded experimentally.) However, in the future an experiment can be done at the 
conditions of very strong screening (i.e. at high salt concentration) with all-purine or all- 
pyrimidine ss-DNA sequences to avoid these contributions. 

The effects of the protein-induced kinks on the ds-DNA elasticity can be studied by 
performing the DNA-stretching experiment at various concentrations of the DNA-binding 
proteins. 

The authors acknowledge valuable discussions with J.-C. Meiners, S. Blumberg, Y. Rabin, 
J. Kierfeld, and T.A. Witten. 
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